Measurement of the cosmic ray antiproton/proton flux ratio at TeV energies with the 

ARGO-YBJ detector. 
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Cosmic ray antiprotons provide an important probe to study the cosmic ray propagation in the 
interstellar space and to investigate the existence of dark matter. Acting the Earth-Moon system 
as a magnetic spectrometer, paths of primary antiprotons are deflected in the opposite sense with 
respect to those of the protons in their way to the Earth. This effect allows, in principle, the search 
for antiparticles in the direction opposite to the observed deficit of cosmic rays due to the Moon 
(the so-called 'Moon shadow'). 

The ARGO-YBJ experiment, located at the Yangbajing Cosmic Ray Laboratory (Tibet, P.R. 
China, 4300 m a.s.L, 606 g/cm 2 ), is particularly effective in measuring the cosmic ray antimatter 
content via the observation of the cosmic rays shadowing effect due to: (1) good angular resolution, 
pointing accuracy and long-term stability; (2) low energy threshold; (3) real sensitivity to the 
geomagnetic field. 
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Based on all the data recorded during the period from July 2006 through November 2009 and on 
a full Monte Carlo simulation, we searched for the existence of the shadow cast by antiprotons in 
the TeV energy region. No evidence of the existence of antiprotons is found in this energy region. 
Upper limits to the p/p flux ratio are set to 5% at a median energy of 1.4 TeV and 6% at 5 TeV 
with a confidence level of 90%. In the TeV energy range these limits are the lowest available. 

PACS numbers: 14.20. Dh;13.85.Tp;96.50.S-;96.50.sd;95.85.Ry 



I. INTRODUCTION 
A. Cosmic ray antiproton production 

Very high energy cosmic ray (VHE CR) antiprotons 
are an essential diagnostic tool to approach the solution 
of several important questions of cosmology, astrophysics 
and particle physics, besides studying fundamental prop- 
erties of the CR sources and propagation medium. The 
enigma of the matter /antimatter asymmetry in the lo- 
cal Universe, that of the existence of antimatter regions, 
the search for signatures of physics beyond the standard 
model of particles and fields, as well as the determina- 
tion of the essential features of CR propagation in the 
insterstellar medium, these are only a few research top- 
ics that would greatly benefit from the detection of VHE 
antiprotons (see for example [HIH)- 

First of all, the observation of p abundance in the CR 
flux is a key to understand CR propagation. In fact, an- 
tiprotons are produced by standard nuclear interactions 
of CR nuclei with the interstellar medium, the informa- 
tion coming from these spallation processes being com- 
plementary to that achievable from secondary nuclei like 
Li, Be and B or secondaries of iron. It should be noticed 
that antiprotons mostly trace the propagation history of 
protons pp— > p ppp, unlike the other spallation products, 
which may come from heavier nuclei. Also secondary p 
represent a background flux that must be carefully de- 
termined to take out a primary p component due to any 
hypothetical exotic signal. 

The observed amount of antiprotons in CRs is still far 
from being figured out. Detailed calculations show that 
there is no model capable of accurately describing alto- 
gether B/C and sub-Fe/Fe ratios, spectra of p, He, p, e + , 
e - and diffuse 7-rays. In fact, conventional models with- 
out reacceleration fail in reproducing both B/C ratio and 
p flux at the same time Q ■ Diffusive reacceleration mod- 
els naturally reproduce secondary/primary nuclei ratio in 
CRs but produce too few antiprotons jfj|. The introduc- 
tion of a break in the diffusion coefficient @ would lead 
to consistent results, but it is not theoretically justified 
so far, still resulting as an ad hoc assumption [4]. Some 
models taking into account Galactic convective wind and 
stochastic reacceleration may reproduce both antiproton 
flux and secondary /primary ratio @. 
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The estimates of the secondary p flux are suffering from 
uncertainties on models and parameters of particle prop- 
agation in the Galaxy, CR spectrum and composition, 
details of the nuclear cross sections for p production, an- 
nihilation and scattering and, finally, on the heliosphcric 
modulation [j| ||. On the contrary, there is a general 
agreement in the calculation of the secondary high-energy 
p flux falloff: around 50 GeV the intensity decreases by 
about 3 orders of magnitude below the maximum. 

Recent measurements of the antiproton flux up to 
about 180 GeV by the PAMELA satellite 0,0] are consis- 
tent with the conventional CR model, in which antipro- 
tons are secondary particles yielded by the spallation of 
CR nuclei over the interstellar medium. Nevertheless, 
given the current uncertainties on propagation parame- 
ters, exotic models of primary p production cannot be 
ruled out [8, 9]. As an example, recent calculations sug- 
gest that the overall PAMELA p and e+ data UBS and 
Fermi e + + e~ data [ll| can be reproduced taking into 
account a heavy dark matter particle (M >10 TeV) that 
annihilates into W + W~ or hh [H|. This scenario implies 
that the p/p ratio, consistent with the background of sec- 
ondary production up to about 50 GeV, increases rapidly 
reaching the 10~ 2 level at about 2 TeV. CR antiprotons, 
as well as positrons, are therefore considered as prime 
targets for indirect detection of dark matter [13l - [l5j . 

But it has been also suggested that the PAMELA 
positron data may be a natural consequence of the stan- 
dard scenario for the origin of galactic CRs, if secondary 
c + (and e _ ) production takes place in the same region 
where CRs are being accelerated [llj]. Since this is a 
hadronic mechanism, an associated rise of the p/p ra- 
tio is predicted at energies > 1 TeV [13]. Therefore, the 
high-energy range of the antiproton spectrum may reveal 
important constraints on the physics of the CR acceler- 
ation sites. 

Antiprotons can be produced from primordial black 
holes evaporation [lj| [l9| or in antigalaxies (2p| - [23l | . In 
particular, it seems possible that mechanisms exist that 
could produce the formation of separated antimatter do- 
mains during the cosmic evolution, thus allowing the vis- 
ible Universe to be globally matter-antimatter symmetric 
[241 ]. In addition, the possibility exists to have antimat- 
ter confined into condensed bodies like antistars in our 
Galaxy [25Ll26|. 

In theories in which matter and antimatter are present 
in equal amounts in spatially separated domains of sur- 
vivable size it is expected that the p/p ratio should 
increase with energy in the framework of the energy- 
dependent confinement model for CRs in the Galaxy. In 
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a simple "leaky box" model the energy spectrum is solely 
determined by the balance between generation at the 
source and escape from the Galaxy. If the source spec- 
trum is proportional to E~", the equilibrium spectrum 
of CRs inside the source region (i.e. inside the Galaxy) 
would be cx E~( a+S \ due to the energy- dependent leak- 
age of the source. Indeed, according to recent measure- 
ments of B/C ratio in the primary CRs up to about 50 
GeV/nucleon, the residence time of CRs in the Galaxy 
can be described by a power law in energy or rigidity 
oc R~ 5 , where S ~0.6 [5]. As a result, the energy spec- 
trum of the CRs leaked out from the antigalaxy has a 
spectrum cx E~ a . If we assume that there exists a general 
acceleration mechanism for generating CRs which acts in 
both galactic and extragalactic sources to give an univer- 
sal source spectrum E _a , the extragalactic CR spectrum 
should reflect it. Thus, if antiprotons are assumed to be 
both primary and extragalactic, we should observe the 
source spectrum of antigalaxies oc E~ a and the expected 
p/p ratio should increase with energy: p/p oc E s [2ll ? 
, [22| ■ As a consequence, the antiproton fraction could in- 
crease up to about 1% around 500 GeV or even to 50% in 
the multi-TeV energy range with important observational 
implications, being, at these energies, the background of 
secondary antiprotons well below this prediction. 

At high energies (> 100 GeV) the main observable re- 
lated to the residence time of CRs in the Galaxy is the 
large-scale anisotropy in their arrival direction that is 
known to be strictly related to the diffusion coefficient. 
Measurements give the amplitude of the first angular har- 
monic of anisotropy at the order 10 -3 in the energy range 
10 11 to 10 14 eV where the most reliable data are avail- 
able [13, HI]- The data on CR anisotropy are consistent, 
within a factor of about 3, with a diffusion coefficient 
increasing with energy oc E 3 , as predicted in models 
including stochastic reacceleration by Kolmogorov-type 
hydromagnetic turbolence (2nd order Fermi acceleration) 
Q. The measurement of the p/p ratio at high energies 
may be useful to constrain models for p production and 
for the confinement of CRs, even if it is not straightfor- 
ward to infer the propagation parameters, as the diffusion 
index 5, since they are partially degenerate with source 
parameters [l7j . 

The first antiproton upper limits in the high energy 
region have been obtained by Stephens in 1985 exploiting 
the observed charge ratio of muons at sea level. The 
presence of p dilutes the charge ratio The limits 

thus derived for the p/p ratio are 7%, 17%, 10% and 14% 
respectively for the energy intervals 0.1 - 0.2 TeV, 1.0 - 
1.5 TeV, 10 - 15 TeV and >30 TeV Q. 

In addition, deeper measurements of the p/p ratio at 
high energies have been performed exploiting the Earth- 
Moon system as a magnetic spectrometer able to disen- 
tangle, in principle, the deflection of protons from that 
of antiprotons in the geomagnetic field. 



B. The Moon shadowing effect 

Since the Moon has an angular radius of about 0.26°, 
it must cast a shadow in the nearly-isotropic CR flux 
(the so-called shadow of the Moon). As first suggested 
by Clark in 1957 [Hj], the shadowing of CRs from the 
direction of the Moon is useful in measuring the angular 
resolution of an air shower array directly, without need- 
ing Monte Carlo (MC) simulations. In fact, the shape 
of the shadow provides a measurement of the detector 
point spread function, and its position allows the check 
of possible pointing biases. 

In addition, due to the geomagnetic field (GMF), posi- 
tively charged particles are deflected by an angle depend- 
ing on the primary CR energy [3(| ■ This effect produces 
a displacement of the shadow towards the West with re- 
spect to the Moon position and smears the shape in the 
East- West direction, especially at low energies. The ob- 
servation of the displacement of the Moon provides a di- 
rect calibration of the relation between shower size and 
primary energy [30| . 

The same shadowing effect can be seen in the direc- 
tion of the Sun. Nevertheless, the interpretation of the 
shadow phenomenology is more complex. In fact, the 
displacement of the shadow from the apparent position 
of the Sun could be explained by the joint effects of 
the GMF and of the Solar and Interplanetary Magnetic 
Fields (SMF and IMF, respectively) whose configuration 
considerably changes with the phases of the solar activity 
cycle (HI, M. 

Linsley [33| and Lloyd-Evans [3J] in 1985, following 
a Watson's suggestion, independently explored the pos- 
sibility to use the Moon or Sun shadows as mass spec- 
trometers in order to measure the charge composition of 
CR spectrum. In particular Linsley first discussed the 
idea to measure the CR antiprotons abundance exploit- 
ing the separation of the proton and antiproton shadows. 
In 1990 Urban et al. [35j carried out detailed calculation 
of this effect proposing this method as a way to search 
for antimatter in primary CR at the TeV energies. 

The GMF should deflect the antimatter component in 
the CRs in opposite direction with respect to the mat- 
ter component. Therefore, if protons are deflected by 
the GMF towards East, antiprotons are deflected to- 
wards West. If the energy is low enough and the angular 
resolution is adequate we can distinguish, in principle, 
two shadows, one shifted towards West due to the pro- 
tons and the other shifted towards East due to the an- 
tiprotons. At high energy (> 10 TeV) the magnetic de- 
flection is too small compared to the angular resolution 
and the shadows cannot be disentangled. At low energy 
(sslOO GeV) the well deflected shadows are washed out 
by the poor angular resolution, thus limiting the sensi- 
tivity. Therefore, there is an optimal energy window for 
the measurement of the antiproton abundance. 

In 1991, the CYGNUS collaboration HI first observed 
the CR shadowing effect measuring a deficit of 4.9 stan- 
dard deviations (s.d.) in the CR background by super- 
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posing the Moon and Sun data at an energy of about 50 
TeV. In the same year also the EAS-TOP experiment ob- 
served the shadowing effect due to the Moon and the Sun 
on the 100 TeV CRs flux with a significance of about 2.7 
s.d. (37| . In the following years this effect has been con- 
firmed by other EAS-arrays (CASA-MIA0I, HEGRA 
[H, GRAPES gjj). 

The first observations of a shadowing effect had to wait 
for the results of the CYGNUS and EAS-TOP experi- 
ments in 1991. There are mainly two reasons for this long 
delay, first the poor angular resolution of EAS-arrays in 
comparison with the angular radius of the Moon or Sun 
(^0.26°). Indeed, only at the beginning of the 90s the an- 
gular resolutions of EAS-arrays reached the 1 deg level. 
Second, due to the high energy threshold (~100 TeV) 
of the experiments the statistical significance of the ob- 
servations was small and the position of the shadow not 
affected by the GMF. Therefore, the deficit of the count- 
ing rate was observed as a function of the angular dis- 
tance from the Moon position, without any information 
on the East- West asymmetry and consequently, without 
any possibility to study the CR antimatter content. 

In 1993, the Tibet AS7 experiment measured both the 
Moon and Sun shadows with an energy threshold low 
enough (about 10 TeV) to allow a 2-dimensional study of 
the effect. In particular, they observed for the first time 
a westward displacement of the Moon shadow (0.16° at 
the 7.1 a level) from its actual position. With this result 
the first upper limit to the p/p ratio with this technique 
was set at about 30% [4lJ. Afterwards, the collaboration 
set an upper limit to the p/p at 10 TeV at about 10% 
42]. 

The CR shadowing effect has been observed also in 
the high energy muon distribution with underground de- 
tectors (SOUDAN-2 El, MACRO Hlli|, L3+C Q, 
BUST |47|, MINOS jH). Recently, the ICECUBE ex- 
periment observed the Moon shadow in the Southern 
hemisphere with a statistical significance of more than 
10 s.d. (49|. Upper limits to the p/p ratio have been set 
by the MACRO (48% at 68% c.l. at about 20 TeV) [H 
and L3+C (11% at 90% c.l. around 1 TeV) [H collabo- 
rations. 

High sensitivity observations of the Moon shadow have 
been recently reported in the multi-TeV energy region by 
an upgraded version of the Tibet AS7 array [50| and by 
the MILAGRO collaboration [IJ. While the Tibet AS7 
experiment set a limit to the p/p ratio to 7% at 90% c.l. 
at about 3 TeV, so far the MILAGRO collaboration did 
not publish any result on the antimatter search with the 
Moon shadow. 

The ARGO-YBJ experiment is particularly effective in 
measuring the CR antimatter content via the observation 
of the CRs shadowing effect due to: (1) good angular 
resolution, pointing accuracy and long-term stability; (2) 
low energy threshold; (3) real sensitivity to the GMF 
due to the absence of any systematic shift in the East- 
West direction. In this paper we report the measurement 
of the p/p ratio in the TeV energy region with all the 



data recorded during the period from July 2006 through 
November 2009. 

The paper is organized as follows. In the Section II the 
ARGO-YBJ detector is described. The description of a 
detailed MC simulation of the Earth-Moon spectrome- 
ter system developed in order to evaluate the deficit of 
events and to calibrate the detector is briefly sketched out 
in Section III. In Section IV the data analysis is outlined 
and the detector performance summarized. The p/p ra- 
tio calculation method is also described in Section IV. C. 
Finally, the results of the data analysis are presented and 
discussed in Section V. A summary of the obtained re- 
sults is given in Section VI. 



II. THE ARGO-YBJ EXPERIMENT 
A. The detector 

The ARGO-YBJ experiment, located at the YangBa- 
Jing Cosmic Ray Laboratory (Tibet, P.R. China, 4300 
m a.s.L, 606 g/cm 2 ), is constituted by a central carpet 
~74x 78 m 2 , made of a single layer of Resistive Plate 
Chambers (RPCs) with ^93% of active area, enclosed 
by a guard ring partially instrumented (~20%) up to 
^100x110 m 2 . The RPC is a gaseous detector work- 
ing with uniform electric field generated by two parallel 
electrode plates of high bulk resistivity (10 n f2 cm). The 
intense field of 3.6 kV/mm at 0.6 atm pressure provides 
very good time resolution (1.8 ns) and the high electrode 
resistivity limits the area interested by the electrical dis- 
charge to few mm 2 . The apparatus has a modular struc- 
ture, the basic data acquisition element being a cluster 
(5.7x7.6 m 2 ), made of 12 RPCs (2.85x1.23 m 2 each). 
Each chamber is read by 80 external strips of 6.75x61.8 
cm 2 (the spatial pixel), logically organized in 10 inde- 
pendent pads of 55.6x61.8 cm 2 which represent the time 
pixel of the detector (52[- The read-out of 18360 pads 
and 146880 strips are the experimental output of the de- 
tector. The RPCs are operated in streamer mode by 
using a gas mixture (Ar 15%, Isobutane 10%, TetraFluo- 
roEthane 75%) for high altitude operation [53]. The high 
voltage settled at 7.2 kV ensures an overall efficiency of 
about 96% [HI . The central carpet contains 130 clusters 
(hereafter ARGO-130) and the full detector is composed 
of 153 clusters for a total active surface of ~6700 m 2 . 
The total instrumented area is ^11000 m 2 . 

A simple, yet powerful, electronic logic has been imple- 
mented to build an inclusive trigger. This logic is based 
on a time correlation between the pad signals depending 
on their relative distance. In this way, all the shower 
events giving a number of fired pads N pa( i > N tr i ff in 
the central carpet in a time window of 420 ns generate 
the trigger. This trigger can work with high efficiency 
down to N trig = 20, keeping negligible the rate of ran- 
dom coincidences. The timing calibrations of the pads is 
performed according to the method reported in [55L [56[ . 

The whole system, in smooth data taking since July 
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2006 with ARGO-130, is in stable data taking with the 
full apparatus of 153 clusters since November 2007 with 
the trigger condition N tr j S = 20 and a duty cycle >85%. 
The trigger rate is ^3.5 kHz with a dead time of 4%. 

Once the coincidence of the secondary particles has 
been recorded, the main parameters of the detected 
shower are reconstructed following the procedure de- 
scribed in [30]. In short, the reconstruction is split into 
the following steps. Firstly, the shower core position is 
derived with the Maximum Likelihood method from the 
lateral density distribution of the secondary particles. In 
the second step, given the core position, the shower axis 
is reconstructed by means of an iterative un-weighted 
planar fit able to reject the time values belonging to the 
non-gaussian tails of the arrival time distribution. Fi- 
nally, a conical correction is applied to the surviving hits 
in order to improve the angular resolution. Unlike the 
information on the plane surface, the conical correction 
is obtained via a weighted fit which lowers the contribu- 
tion from delayed secondary particles, not belonging to 
the shower front. 

The analysis reported in this paper refers to events se- 
lected according to the following criteria: (1) more than 
25 strips N strip should be fired on the ARGO-130 car- 
pet; (2) the zenith angle of the shower arrival direction 
should be less than 50°; (3) the reconstructed core posi- 
tion should be inside an area 150 x 150 m 2 centered on the 
detector. After these selections the number of events an- 
alyzed is about 2.5 x 10 11 (about 10 9 inside a 10° x 10° an- 
gular region centered on the Moon position). According 
to simulations, the median energy of the selected protons 
is E 50 wl.8 TeV (mode energy «0.7 TeV). 



III. MONTE CARLO SIMULATION 

The deficit of events in the Moon direction, the ab- 
solute energy calibration and the angular resolution, as 
well as the systematic pointing biases, have been stud- 
ied by comparing the observed Moon shadow charac- 
teristics (East- West and North-South displacements and 
shape) with the expectations from a detailed MC simu- 
lation of the CR propagation in the Earth-Moon system 

With this simulation we also estimated the expected 
antiprotons flux in the opposite CR Moon shadow side, 
as described in Section IV. C. 

From the MC simulation strategy viewpoint, the Moon 
shadow has been treated like an extensive excess signal, 
instead of a lack in the isotropic CRs flux. In other words, 
our simulation deals with the Moon as if it was the source 
of the CRs which it intercepts in reality. 

The simulation has been realized on the basis of the 
real data acquisition time. The Moon position has been 
computed at fixed times, starting from July 2006 up to 
November 2009. Such instants are distant 30 seconds 
each other. For each time, after checking the data acqui- 
sition was effectively running and the Moon was in the 



field of view, extensive primaries are generated with ar- 
rival direction sampled within the Moon disc. For each 
chemical species (p, He, CNO group, Mg-Si group and 
Fe), the number of primaries to be generated is computed 
on the basis of the effective exposure time, according to 
the energy spectrum resulting from a global fit of the 
main experimental data (59|. 

Once the number of CRs expected to be hampered 
by the Moon has been calculated, the charge sign of ev- 
ery primary is inverted and it is propagated back to the 
Moon, the magnetic field bending its trajectory. The 
propagation stops anyway at the Moon distance, giving 
a good approximation of the deflection undergone by the 
CR before reaching the atmosphere. In fact, if we firstly 
consider a positively charged CR arriving to the Earth 
atmosphere and then we invert its charge and its momen- 
tum and threw it back to the space, the two trajectories 
do not overlap because of the numerical approximation. 
Nonetheless, as long as the primary energy is above sev- 
eral tenth of GeV, both trajectories give similar devia- 
tions, the difference ranging from 15% at 50 GeV down 
to 3% above 1 TeV. Further details and results from this 
simulation can be found in [3(3, H3] • 

After accounting for the arrival direction correction 
ought to the magnetic bending effect, the air showers 
development in the atmosphere has been generated with 
the CORSIKA v. 6.500 code The electromagnetic 
interactions are described by the EGS4 package while the 
hadronic interactions above 80 GeV are reproduced by 
the QGSJET-II.03 and the SYBILL models. The low en- 
ergy hadronic interactions are described by the FLUKA 
package. CR spectra have been simulated in the energy 
range from 10 GeV to 1 PeV following the relative nor- 
malization given in 59]. About 10 8 showers have been 
sampled in the zenith angle interval 0-60 degrees. The 
secondary particles have been propagated down to a cut- 
off energies of 1 MeV (electromagnetic component) and 
100 MeV (muons and hadrons). The experimental condi- 
tions (trigger logic, time resolution, electronic noises, re- 
lation between strip and pad multiplicity, etc.) have been 
taken into account via a GEANT4-based code [6l| . The 
core positions have been randomly sampled in an energy- 
dependent area large up to 2-10 3 x 2-10 3 m 2 , centered 
on the detector. Simulated events have been generated 
in the same format used for the experimental data and 
analyzed with the same reconstruction code. 



IV. DATA ANALYSIS 

For the analysis of the shadowing effect, the signal is 
collected within a 10°xl0° sky region centered on the 
Moon position. We used celestial coordinates (right as- 
cension and declination, R.A. and DEC. hereafter) to 
build the event and background sky maps, with 0.1° x0.1° 
bin size. Finally, after a smoothing procedure, the signif- 
icance map, used to estimate the statistical significance 
of the observation, is obtained. 
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FIG. 1: Plot (a): showers firing N s tri P >100 collected around 
the Moon position. The coordinates are R.A. a and DEC. 5 
centered on the Moon position (a m , S m ). The plot (b) shows 
the map projections along the R.A. direction. 

As it can be appreciated in Fig. [TJ the Moon shadow 
turns out to be a lack in the smooth CR signal, observed 
by ARGO-YBJ even without subtracting the background 
contribution nor smoothing the signal. 



A. Moon shadow analysis 

Cosmic rays blocked by the Moon must be as much 
as the background events lying within a region as large 
as the Moon disc. A suitable background estimation is 
therefore a crucial point of the analysis. As it can be 
seen also from the projection along the R.A. from Fig. 
[TJ the background events are not uniformly distributed 
around the Moon, because of the non-uniform exposure 
of the map bins to the CR radiation. The background 
has been estimated with the equi-zenith angle method, 
as described in detail in [30j . 

A significance map of the Moon region is shown in 
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FIG. 2: Significance map of the Moon region observed with 
all events detected by ARGO-YBJ. The event multiplicity is 
25< N atrip < 40 and zenith angle 9 < 50°. The coordinates 
are R.A. a and DEC. S centered on the Moon position (a m , 
S m )- The color scale gives the statistical significance in terms 
of standard deviations. 

Fig. [21 It contains all events belonging to the lowest 
multiplicity bin investigated (25< N s t r i P < 40), collected 
by ARGO-YBJ during the period July 2006 - November 
2009 (about 3200 hours on-source in total). N s t r i P is the 
number of fired strips on the central carpet ARGO-130. 
The significance of the maximum is about 22 s.d.. The 
observed westward displacement of the Moon shadow by 
about 1.5° allows to appreciate the sensitivity of the 
ARGO-YBJ experiment to the GMF. This means that 
a potential antiproton signal is expected eastward within 
1.5° from the actual Moon position (i.e., within 3° from 
the observed Moon position). The median energy of se- 
lected events is E 50 «750 GeV (mode energy « 550 GeV) 
for proton-induced showers. The corresponding angular 
resolution is ~1.6°. 

The large displacement of the shadow is only one el- 
ement of this analysis, the other one being the angu- 
lar resolution which is not adequate in this multiplicity 
range. Indeed, as can be seen from the Fig. [2j the mat- 
ter shadow is visible on the antimatter side with a sig- 
nificance of about 10 s.d., thus limiting the sensitivity to 
the antiproton abundance measurement. We note that 
this is the first time that an EAS array is observing the 
Moon shadow cast by sub-TeV primary CRs. 

B. Detector performance 

The performance of the detector and its operation sta- 
bility have been studied in detail in [30| exploiting the 



CR Moon shadowing effect with all data since July 2006. 

The measured angular resolution is better than 0.5° for 
CR-induced showers with energies E > 5 TeV, in good 
agreement with MC expectations. The Point Spread 
Function of the detector, studied in the North-South 
projection not affected by the GMF, is Gaussian for 
N strip >200, while for lower multiplicities is better de- 
scribed for both MC and data with a linear combination 
of two Gaussian functions. The second Gaussian con- 
tributes for about 20%. 

The long-term stability of the ARGO-YBJ experiment 
has been checked by monitoring both the position of the 
Moon shadow, separately along R.A. and DEC. projec- 
tions, and the amount of shadow deficit events in the pe- 
riod November 2007 - November 2010, for each sidereal 
month and for events with N str i P > 100. As shown in Fig. 
17 of ref. [30(, the position of the Moon shadow turned 
out to be stable at a level of 0.1° and the angular resolu- 
tion stable at a level of 10%, on a monthly basis. These 
results make us confident about the detector stability in 
the long-term observation of the Northern sky. A sys- 
tematic uncertainty of (0.19±0.02)° towards the North 
in the absolute pointing accuracy is observed. The most 
important contribution to the systematics is likely due 
to a residual effect not completely corrected by the time 
calibration procedure. Further studies are under way. 

We have estimated the primary energy of the detected 
showers by measuring the westward displacement as a 
function of the shower multiplicity, thus calibrating the 
relation between shower size and CR energy. The system- 
atic uncertainty in the absolute rigidity scale is evaluated 
to be less than 13% in the range from 1 to 30 TeV/Z, 
mainly due to the statistical one [30]. 



C. The CR p/p flux ratio calculation 

All chemical species of CR, each with its own spectrum, 
contribute to form the Moon shadow signal. Hence, the 
chance of unfolding all contributions relies on MC simu- 
lations, as well as the search for antiprotons demands to 
properly reproduce the Moon shadow signal. The shape 
of the Moon shadow is tightly connected to the primaries 
energy spectrum, which mostly determines the tails of 
the signal. At first, we assume the energy spectrum of 
antiprotons follows a power law dN/dE = k ■ E r , where 
the spectral index 7 is taken to be as large as that of pro- 
tons. An investigation of the dependence on the spectral 
index will follow in the Section V. 

In order to evaluate the CR p/p flux ratio, only the 
projection along the R.A. is important, as it has been 
shown that at Yangbajing the GMF has a non-null effect 
on the CR trajectories only along such direction [301 ]. 
The R.A. projection results from an integration along 
the DEC. direction. 

The deficit counts observed around the Moon projected 
on the R.A. axis (that is on the East- West axis) are shown 
in Fig. |3]for two multiplicity bins compared to MC expec- 
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FIG. 3: Deficit counts measured around the Moon projected 
along the East- West axis for two different multiplicity bins 
(black circles) compared to MC expectations (red squares). 
Events contained in an angular band parallel to the East- West 
axis and centered on the observed Moon position, propor- 
tional to the multiplicity-dependent angular resolution, are 
used (see text). 



tations: 40<N strip <100 and N strip >100, in the upper 
and lower panels, respectively. The vertical axis reports 
the events contained in an angular band parallel to the 
East- West axis and centered on the observed Moon posi- 
tion. The widths of these bands are chosen on the basis 
of the MC simulation so that the shadow deficit is max- 
imized. They turn out to be proportional to the N str j p - 
dependent angular resolution. The widths of these bands 
are ±2.80° and ±1.90°, respectively. 

The data are in good agreement with the MC sim- 
ulation and the observed shadows are shifted westward 
of (-0.75±0.05)° and (-0.30±0.05)°, as expected. A de- 
tailed analysis of this sort of projections as a function of 
the shower size is given in [30] ]. 

The GMF shifts westward the dip of the signal from 
positively charged primaries. Searching antiprotons 
means looking for excesses in the eastern part of the R.A. 
projection, i.e. trying and fitting the Moon shape ex- 
pected from combining CR and antimatter to the shape 
obtained from experimental data. Of course, whichever 
matter-antiprotons combination is obtained, the total 
amount of triggered events must not be changed, so that 
the fitting procedure consists in transferring MC events 
from the CR to the antiprotons shadow and comparing 
the result with data. 

To make such a comparison, we firstly adopted the fol- 
lowing method. We obtained two kinds of Moon shadow, 
cast by all CRs and protons, respectively. After project- 
ing them along the R.A. direction, we used a superposi- 
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tion of several Gaussian functions to describe the deficit 
event distribution in each shadow [5(3]. Four Gaussian 
functions were found to be adequate for fitting both dis- 
tributions within 5° from the Moon disc center. Let us 
name 9 the angular distance from the Moon disc center 
and fm(&) the Gaussian function superposition describ- 
ing the CR shadow. Let J- p {6) be the proton shadow, 
obtained by imposing a given power law spectrum. The 
observed Moon shadow should be expressed by the fol- 
lowing function: 

!moon{0) = (1 - r) f m (9) +rJ- ¥ (0) 
= (l-r)f m (6) + rF p (-6) 

(0 < r < 1) where the first term represents the deficit 
in CRs and the second term represents the deficit in an- 
tiprotons. This function must be fitted to the data to 
obtain the best value of r. 

We also applied a second method to determine the an- 
tiproton content in the cosmic radiation. Without intro- 
ducing functions to parameterize the expectations, we 
directly compared the MC signal with the data. We per- 
formed a Maximum Likelihood fit using the p content as 
a free parameter with the following procedure: 

1. the Moon shadow R.A. projection has been drawn 
both for data and MC. 

2. the MC Moon shadow has been split into a "mat- 
ter" part plus an "antiproton" part, again so that 
the total amount of triggered events remains un- 
changed: 

$Mc(mat) — > $Mc(r;mat + p) 
= (1 - r)<P MC (mat) + <$>mc{P) 

3. for each antiproton to matter ratio, the expected 
Moon shadow R.A. projection &Mc( r ', ma t + P) is 
compared with the experimental one via the calcu- 
lation of the likelihood function: 

B 

C{r) =Y,Niln[Ei(r)] - E t (r) - ln(JV,!) 
i=l 

where Nj is the number of experimental events in- 
cluded within the i-th bin, while Ei(r) is the num- 
ber of events expected within the same bin, which 
is calculated by adding the contribution expected 
from MC ($Mc( r ; mat + p)) to the measured back- 
ground. 

Both methods described above give results consistent 
within 10%. 



V. RESULTS AND DISCUSSION 

The optimal energy windows for the measurement 
of the antiproton abundance in CRs is identified by 



the following multiplicity ranges: 40<N str i p <100 and 
N strip >100. In the former bin the statistical significance 
of Moon shadow observation is 34 s.d., the measured an- 
gular resolution is the proton median energy is 1.4 
TeV and the number of deficit events about 183000. The 
shadow is shifted of about 1° westward. In the latter 
multiplicity bin the significance is 55 s.d., the measured 
angular resolution ^0.6°, the proton median energy is 5 
TeV and the number of deficit events about 46500. The 
shadow is shifted of about 0.4° westward. 

There exists, however, no evidence indicating deficits 
of CRs at the opposite positions around 0—1° and 0.4° 
in the eastward direction, corresponding to the particles 
with negative charge (anti-matter) such as p, He, C,... 
and Fe, if any. 

Therefore, we applied both methods described in the 
previous section to evaluate upper limits to the CR p/p 
flux ratio. The r parameter which best fits the expecta- 
tions to the data turns out to be always negative, i.e. it 
assumes non-physical values throughout the whole energy 
range investigated. With a direct comparison of the R.A. 
projections, the r-values which maximize the likelihood 
are: -0.076±0.040 and -0.144±0.085 for 40<N st „ p <100 
and N s trip >100, respectively. The corresponding upper 
limits with 90% confidence level (c.L). according to the 
unified Feldman & Cousins approach [62|, are 0.034 and 
0.041, respectively. 

Since the anti-shadow was assumed to be the mirror 
image of the proton shadow, we assume for the antipro- 
tons the same median energy. The p/p ratio is Q(p) 
= l/f p - <&(p)/$ (matter), therefore, being the assumed 
proton fraction f p = 73% for 40<N str i p <100 and f p = 
71% for N strip >100 [H^l, we are able to set the follow- 
ing upper limits at 90% c.L: 0.05 for 40<N st „ p <100 
and 0.06 for N s trip >100. Notice that the two values are 
similar, in spite of the different multiplicity interval. It 
is a consequence of the combination of the two opposite 
effects of the angular resolution and of the geomagnetic 
deviation. 

In Fig. 2] the ARGO-YB J results are compared with 
all the available measurements. The energy bin is 34% 
around the median energy for each multiplicity interval. 
The solid curves refer to a theoretical calculations for a 
pure secondary production of antiprotons during the CR 
propagation in the Galaxy by Donato et al. [14]. The 
curves were obtained using the appropriate solar mod- 
ulation parameter for the PAMELA data taking period 
[8j. The long-dashed line refer to a model of extragalac- 
tic primary p production [20I, |23| . The rigidity-dependent 
confinement of CRs in the Galaxy is assumed oc R -<5 , be- 
ing R the rigidity, and 5 = 0.6. We note, however, that 
this curve has been normalized by authors to the low 
energy p/p measurements carried out in 1980s. Recent 
measurements show that the sub- TeV p/p flux ratio is 
about a factor 10 lower. The dotted line refers to a pos- 
sible contribution of antiprotons from the annihilation of 
a heavy dark matter particle T2|. The short-dashed line 
shows the calculation by Blasi and Serpico 17J for sec- 
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FIG. 4: The antiproton to proton flux ratio obtained with the ARGO-YBJ experiment compared with all the available mea- 
surements. The solid curve refers to a theoretical calculations for a pure secondary production of antiprotons during the CR 
propagation in the Galaxy by Donato et al. [l4| . The long-dashed line refer to a model of extragalactic primary p production 
[201. [231]. The rigidity-dependent confinement of CRs in the Galaxy is assumed to be tx K~ s , where 8 — 0.6. The dotted line 
refers to the contribution of antiprotons from the annihilation of a heavy dark matter particle 12]. The short-dashed line 
shows the calculation by Blasi and Serpico [171 ] for secondary antiprotons including an additional p component produced and 
accelerated at CR sources. 



ondary antiprotons including an additional p component 
produced and accelerated at CR sources. 

Two sources of systematic errors have been investi- 
gated: the presence of electrons in the cosmic radiation 
and the unknown antiproton spectral index. As discussed 
in 30] the absolute rigidity scale uncertainties associated 
to the CR chemical composition and to different hadronic 
models in the MC calculations are about 7% and 12% re- 
spectively. 

Like antiprotons, electrons are supposed to shift west- 
ward, giving a Moon shadow opposite to that produced 
by positively charged CR. Nonetheless, the Moon shadow 
provided by electrons is expected to be slightly different 
from that of antiprotons. Firstly, for a given multiplic- 
ity range the median energy of electron-induced show- 
ers is 30%-40% less than that of CR showers, i.e. the 
shadow dip is expected to be further displaced 30% - 
40%. Then, the angular resolution is better for electron 
primaries than for CRs (~30%). Last, the electron flux 
at TeV energies is less than 10" 3 of CR flux (63|. As 
a consequence, we estimate that the systematic uncer- 
tainty due to misinterpreting electrons as antiprotons is 
below 10%. 

Since the spectral index of antiprotons is unknown, 
there is no reason to assume the proton spectral index. 



Many unknown factors contribute to its value, mostly re- 
lated to the diffusion coefficient inside galaxy. To inves- 
tigate this point, primary antiprotons are assigned differ- 
ent spectral indices, 7= 2.0, 2.2, 2.4, 2.6, 2.8, 3.0. Results 
for the investigated multiplicity intervals are summarized 
in tablelU The limits of the antiproton/proton ratio varies 
of 20%-30% with the spectral index. 



TABLE I: The effect of different spectral indices 



index 


90% upper limit(40-100) 


90% upper limit(>100) 


2.0 


3% 


4% 


2.2 


4% 


4% 


2.4 


4% 


4% 


2.6 


5% 


5% 


2.8 


5% 


7% 


3.0 


6% 


7% 



VI. CONCLUSIONS 

The ARGO-YBJ experiment is observing the Moon 
shadow with high statistical significance at an energy 
threshold of a few hundreds of GeV. Using all data col- 
lected until November 2009, we set the upper limits on 
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the p/p flux ratio to 5% at an energy of 1.4 TeV and 6% 
at 5 TeV with a confidence level of 90%. 

In the few- TeV energy range the ARGO-YBJ results 
provide the strongest p/p limits obtained to date, useful 
to constrain any primary antiproton production model 
which foresees high fluxes at TeV energies. As dis- 
cussed in Section IV. A the main limiting factor in the 
p/p ratio measurement exploiting the Moon shadow tech- 
nique is the angular resolution. The new generation of 
EAS-arrays under construction (HAWC, LHAASO) is ex- 
pected to improve the angular resolution by a factor of « 
3 in the TeV energy range IH[ • Taking into account 
also an expected increase of the effective area by at least 
a factor of 3, we expect that in the next future the sen- 
sitivity to the p/p ratio could be lowered by a factor of 
10 in the TeV energy range. 
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